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Abstract. 

Colloidal suspensions and polyclcctrolytc solutions containing multivalent countc- 
rions can exhibit some very counter-intuitive behavior usually associated with the low 
temperature physics. There are two particularly striking phenomena resulting from 
strong electrostatic correlations. One is the like-charge attraction and the second is 
the polyion overcharging. In this contribution we will concentrate on the problem of 
overcharging. In particular we will explore the kinetic limitation to colloidal charge 
inversion in suspensions containing multivalent counterions. 
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1. Introduction 

Colloidal suspensions and polyelectrolyte solutions containing multivalent counterions 
can exhibit some very curious electrostatic behavior ffl. It is found that under some 
circumstances two like-charged polyions inside suspension can actually attract one 
another §, |, |, |, |, 0, |, |, Q EH S E3. S S EH- The counterion mediated 



attraction is responsible for the DNA compaction inside the bacteriophages, viruses that 
infect bacteria flTj], |1| , and for the organization of eukaryotic cytoskeleton . Another 



"strange" electrostatic behavior which can occur in suspensions containing multivalent 
counterions is the reversal of the electrophoretic mobility fl], 20, 21, £22], |23], |25| . The 



first thing that is learned in a course on electrostatics is that the force produced by the 
electric field on a charged particle is 

F = QE . (1) 

Thus, a positively charged particle, Q > 0, is expected to move in the direction 
of the applied field while a negatively charged particle, Q < 0, will move in the 
direction opposite to the field. This simple picture, however, breaks down inside 
colloidal suspensions of low dielectric solvent or even in aqueous suspensions containing 
multivalent counterions. The reason for the violation of the "simple" physics learned 
in high-school are the strong electrostatic many-body interactions between the colloidal 
particles and the counterions. The reversal of electrophoretic mobility can be understood 
as a combination of two electrostatically driven mechanisms. Strong electrostatic 
interaction between colloids and counterions leads to formation of polyion-counterion 



complexes []26L |27|, 28fl . The existence of counterion condensation has been known 



for over thirty years |29|, |30], pl| , the general phenomenon is, however, much older 
than this and can be traced to the pioneering work of Bjerrum on ionic association 
inside electrolyte solutions almost 80 years ago p2| . In aqueous suspensions with only 



monovalent counterions, the net charge of complexes is of the same sign as the bare 
charge of polyions. 

If the solvent is water and the counterions are monovalent, the electrostatic 
interactions between the condensed counterions can be neglected 0, and the 
simplest Poisson-Boltzmann theory is sufficient to describe the polyion-counterion 



complexation |26|, In aqueous suspensions containing multivalent counterions or 
in suspension of low dielectric solvents, the electrostatic energy between the condensed 
counterions is significantly larger than the thermal energy and the electrostatic 
correlations between the condensed counterions can no longer be neglected. These 
electrostatic correlations can lead to colloidal overcharging i.e. the net charge of the 
complex is of opposite sign to the charge of the bare polyion. The overcharged colloid 
will then move in the "wrong" direction with respect to the applied electric field [|TJ. [25 
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2. Overcharging 

To understand the phenomenon of overcharging we shall start by studying a very simple 
model. Consider a sphere of radius a and fixed charge — Zq distributed uniformly over 
its surface. We would like to know how many point-like a-valent counterions, each 
of charge acq, should be placed on top of this sphere in order to minimize the total 
electrostatic free energy |34|, ^j, 24|. When we say "counterions" we have in mind 



both simple multivalent ions such as Ca ++ , as well as more complicated micelle-like 
aggregates with a significantly higher than one. 

The free energy of a complex can be written as 

lea ea 

The first term is the self energy of the charged sphere, the second term is the electrostatic 
energy of interaction between the sphere and n condensed a-ions, and the last term is 
the electrostatic energy of repulsion between the condensed counterions. To calculate 
the free energy of repulsion, it is convenient to express F£ a in terms of the free energy 
of a one component plasma (OCP), n a-ions on the surface of a sphere with a uniform 
neutralizing background, F® CP . The free energy of a spherical OCP can be written as 

P ocp P aa a 2 nV a 2 n 2 q 2 

b n = b n 1 o • \ 6 > 

ea lea 

Substituting Eq. (|3|) into Eq. (0) the electrostatic free energy of a polyion-counterion 
complex becomes, 

_ (Z - anfq 2 OCP 
E n - — + F n . (4) 

In the strong coupling limit, corresponding to multivalent counterions or solvents of low 
dielectric permittivity, the free energy of the OCP is well approximated by the free 
energy of the low temperature phase corresponding to a triangular Wigner crystal, 

F °^-M^. ,5, 

where M is the Madelung constant. For weaker couplings, the expression for the F® CP 
can be obtained from the fits to the Monte Carlo data |35| . For concreteness we shall 
use M = 1.106, the value appropriate for a planar Wigner crystal [|l|]. 

The effective charge of a polyion-counterion complex, in units of — q is 

Z eff = Z -an , (6) 

The optimum number of condensed counterions is determined from the minimization of 



the total electrostatic free energy. We find [22], [24 , 



where 



7= 3m7^ (8) 
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We see that the optimal charge of a polyion-counterion complex is of opposite 
sign to the bare colloidal charge, i.e. the complex is overcharged. Inside the colloidal 
suspension containing multivalent counterions or solvents of low dielectric permittivity 
the electrophoretic mobility can, therefore, be reversed. 

Some care, however, must be taken in extrapolating the results of this simple model 
to real systems. While we have treated the counterions as condensed on top of the 
sphere, this is clearly not the case for real colloidal suspension. Instead the associated 
counterions form a layer around a colloidal particle which can be some nanometers wide. 
The presence of simple electrolyte also strongly affects the net charge of the polyion- 
a-ion complex. Furthermore, the complex formation is a kinetic phenomenon requiring 
a counterion to overcome an energy barrier in order to join the already overcharged 
complex. 



3. The overcharging potential 

In the previous section we found that the minimum of the total electrostatic free 
energy of a polyion-a-ion complex corresponds to an overcharged state. However, for a 
counterion to join an already overcharged complex it must overcome an energy barrier. 
The waiting time for a thermal fluctuation of sufficient strength necessary to drive a 
counterion over an activation barrier scales exponentially with the height of the barrier. 
There is, therefore, a kinetic limitation to the degree of overcharging which can prevent 
a thermodynamically optimum state from being reached on experimental time scale. To 
explore this further we have to construct an effective interaction potential between a 
complex and a counterion separated by distance r. 

The work necessary to bring a counterion from infinity to join a complex containing 
n a-ions is 

W = ^. (9) 

an 

We define the reduced electrostatic potential of a counterion on the surface of the 
complex as <f (a) = (3W , where j3 = l/k B T. Differentiating Eq. (|) we find 

. . (Z — an)\ B a 3Ma 2 \/n , . 

<p(a) = - K - ^ — , 10 

a 4a 

where X B = q 2 /tk B T . The first term of Eq. (0) is the electrostatic energy of interaction 
between a uniform spherical charge and an a-ion, while the second term is due to 
electrostatic correlations between the a-ions. In the strong coupling limit correlational 
contribution to the interaction potential decay exponentially fast with the separation 
from the polyion surface |7], [l^, [36|. The characteristic length is set by the average 
separation between the condensed counterions. More specifically we can approximate 
the reduced interaction potential by 

^ = (Z - an)\ B a _ 3Ma> e _ (r _ a)/e (n) 
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The decay of the correlational contribution is governed by the characteristic length £ 



which in the strong coupling limit is well approximated by [7L 15, Bq 



f = M' (12) 

where G is the reciprocal lattice vector of a triangular Wigner crystal of condensed 



counterions. Due to strong coupling between the condensed counterions, Eq. (|T2|) should 
remain a good approximation even significantly above the crystallization temperature. 
For a triangular Wigner crystal, 

|G| ^ (13) 

where b is the lattice spacing 

b = (14) 

and a = n/Ana 2 is the surface density of condensed counterions. Substituting Eqs. ([13|) 
and fli~4]) into Eq. (|12|) , the decay length is found to be 



{ = wr (15) 

We are now in possession of the electrostatic potential which will allow us to study the 
kinetics of overcharging. 

4. Kinetics of overcharging 

For n < Z/a, the electrostatic potential between a counterion and a complex is purely 
attractive favoring further counterion condensation. Inside an electrolyte solution this 
tendency towards polyion-counterion association is opposed by the loss of entropy 
resulting from the confinement of condensed counterions near the colloidal surface. Here, 
however, we shall not be concerned with the role of entropy |TJ. 

For n > Z/a the interaction potential has two minima, one located at r = a and 
the second at n = oo. For Z/a < n < n* the r = a minimum is the dominant one, 
while for Z > n* the global minimum changes to r = oo. The value of n* corresponds 
to the number of condensed counterions which minimize the electrostatic free energy of 
the complex Eq. (f|), 

Z Z \ r p 

n* = 2Z£ . (1 6 ) 

a 

In the case of trivalent counterions the energy barrier that a counterion needs to 
overcome in order to join a complex which already contains n* condensed a-ions is 
less than 2/cgT, Fig 1. Thus, for trivalent counterions there is no kinetic hindrance to 
reaching the optimum overcharged state. 

We next look at the height of the activation barrier as a function of the counterion 
valence, Fig. 2. It is clear that the height of the activation barrier grows rapidly with the 
increased valence of the a-ions. In particular we see that for a = 10 the activation barrier 
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Figure 1. The reduced interaction potential between a complex of Z = 4000, a = 1000 
A, n = n* condensed trivalent counterions, and a trivalent counterion located at 
distance r from the center of colloid. 




1 6 11 16 21 26 

Valence of counterions 

Figure 2. The height of the activation barrier that an a- ion must overcome to join 
an optimally overcharged complex composed of a colloid with Z = 4000, a — 1000 A 
and n — n* condensed a-ions. 
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is already some 10fc#T which is probably the maximum height that a counterion can 
overcome on a reasonable experimental time scale. Thus, the process of overcharging 
by the ct-ions with a > 10 will be kinetically controlled. For example, from Eq. (0) 
we see that the optimal state of overcharging of a colloidal particle of Z = 4000 and 
radius a = 1000 A by micelles with a = 25 corresponds to Z* e ^ = —271. In practice, 
though, the process of overcharging will come to a stop when the barrier height reaches 
about 10k bT, implying that the complex will stop growing when the net charge is only 

Zeff = -70. 

5. Conclusion 

In this contribution we have explored the kinetic limitation to overcharging. We find 
that kinetics does not play an important role for overcharging by simple multivalent 
counterions, so that the state of optimal overcharging, Eq. (^), is accessible within an 
experimental time scales. On the other hand, we find that the activation barrier grows 
rapidly with the valence of counterions, suggesting that the extent of overcharging by 
micelle-like aggregates is largely kinetically controlled. 

The kinetic limitation to overcharging might also be important for the formation 
of the DNA-cationic lipid complexes. The problem of a reliable and safe mechanism for 
gene delivery is particularly pressing in view of the current medical applications. Strong 
electrostatic repulsion between a DNA and a cellular membrane inhibits transfection of 
a naked DNA into the cell. A way to overcome this difficulty is through the formation of 
overcharged complexes between the DNA and the cationic liposomes |37j , |38| , |39| , |40| , |4l 



These lipoplexes having a net positive charge are attracted to the cellular membrane, 
facilitating the genetic transfection. 

Finally, the presence of a simple electrolyte will have a strong influence on the 
overcharging. It has been demonstrated that for sufficient concentration of a-ions, 
monovalent salt favors overcharging ||2| [l[ . In fact in the presence of simple electrolyte 
the thermodynamic state of optimum overcharging corresponds to the charge inversion 
of as much as 100%. This should be contrasted with the result of Eq. (0), which shows 
that in the absence of salt, the effective charge of a complex scales as a square root of 
the bare charge. The presence of salt will also lower the height of the activation barrier 
reducing the kinetic hindrance to overcharging. 
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